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Abstract. We study the influence of an imperfect structure of a crystalline undulator 
on the spectrum of the undulator radiation. The main attention is paid to the 
undulators in which the periodic bending in the bulk appears as a result of a regular 
(periodic) surface deformations. We demonstrate that this method of preparation of 
a crystalline undulator inevitably leads to a variation of the bending amplitude over 
the crystal thickness and to the presence of the subharmonics with smaller bending 
period. Both of these features noticeably influence the monochromatic pattern of the 
undulator radiation. 



1. Introduction 

In this paper we present the qualitative and quantitative analysis of the influence of 
imperfect structure of a crystalline undulator (in particular, the influence of the variation 
of the bending amplitude over the crystal thickness due to the stress applied to its 
surface) on the spectral distribution of the radiation. The parameters of crystalline 
undulators (including the types and lengths of crystals, the periods of bending, the 
positron energies, and the range of photon energies) as well as the methods of preparation 
of periodically bent structures discussed below in the paper correspond to those which 
are available for the experiments to be carried out within the PECU project [T]. 

A periodically bent crystal together with a bunch of ultra-relativistic charged 
particles which undergo planar channeling constitute a crystalline undulator. In such a 
system there appears, in addition to the well-known channeling radiation, the undulator 
type radiation which is due to the periodic motion of channeling particles which follow 
the bending of the crystallographic planes[2j [3]. The intensity and characteristic 
frequencies of this radiation can be varied by changing the beam energy and the 
parameters of bending. In the cited papers as well as in the subsequent publications 
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(see the review Ref. [I] and the references therein) a feasibility was proven to create a 
short-wave crystalline undulator that will emit intensive monochromatic radiation when 
a pulse of ultra-relativistic positrons is passed through its channels. More recently, it was 
demonstrated [5] that the undulator based on the electron channeling is also feasible. 
A number of corresponding novel numerical results were presented to illustrate the 
developed theory, including, in particular, the calculation of the spectral and angular 
characteristics of the new type of radiation. 

Although the operational principle of a crystalline undulator does not depend on the 
type of a projectile below in this paper we will consider the case of a positron channeling. 
Under certain conditions [2, [3] , an ultra-relativistic positron, which enters the crystal 
at the angle smaller than the Lindhard critical angle [6], will penetrate through the 
crystal following the bendings of its planes. Provided the bending amplitude a greatly 
exceeds the interplanar distance d (see figure [U left) one can disregard the oscillations 
due to the action of the interplanar force, - the channeling oscillations [6]. In this case 
the trajectory of the particle can be associated with the periodic profile of channel 
centerline. The undulator radiation appears as a result of this periodic motion of the 
particle. Thus, the operational principle of a crystalline undulator is the same as for a 
conventional one[7],[8], in which the monochromaticity of the radiation is the result of a 
constructive interference of the photons emitted from similar parts of trajectory. 




Figure 1. Schematic representation of a crystalline undulator with a constant bending 
amplitude a (left panel) and with a varied amplitude a(y) (right panel). Circles 
denote the atoms belonging to neighbouring crystallographic planes (separated by the 
distance d) which are periodically bent with a period A. The dashed curves denote the 
centerlines of the planar channels. 



Usually, when discussing the properties of a crystalline undulator and of its 
radiation, one considers the case of a perfect crystalline undulator. By this term we will 
understand the crystal whose planes are bent periodically following a perfect harmonic 
shape, y(z) = asin(27rz/A), see figured] (left). For clarity, let us stress that we consider 
the case when the quantities d, a and A satisfy the double inequality d <C a <C A. 
Typically, d ~ 1 A, a — 10 . . . 10 2 d and A ~ 10~ 5 . . . 10~ 4 a. The spectral-angular 
distribution from the perfect undulator is characterized by a specific pattern which 
implies that for each value of the emission angle the spectrum consists of a set of narrow, 
well-separated and powerful peaks corresponding to different harmonics of radiation. In 
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principle, the perfect crystalline undulators can be produced by using the technologies 
of growing Sii^Ge^ structures [9]. In this case, by varying the Ge content x one can 
obtain periodically bent crystalline structure [TQl HI]- The technological restrictions 
imposed by this method on the crystalline undulator length is L < 140 . . . 150 /im. 

The periodic bending can also be achieved by making regularly spaced grooves 
on the crystal surface either by using a diamond blade [121 US] or by means of laser- 
ablation. The latter method was used [H] to prepare the Si-based crystalline undulators 
for the PECU experiments. In either case, the regular surface deformation results in 
the periodic pattern of the crystallographic planes bending in the bulk. The question 
which appears in connection with these methods of preparation concerns the quality 
of the periodically bending. Indeed, for a crystal of a finite thickness it is natural to 
expect that the surface deformations, regularly spaced with the period A, result in the 
volume deformations of the same period but of a varied amplitude of bending, a = a(y) 
(see figured] (right)). The latter has maximum value in the surface layer but decreases 
with the penetration distance. Therefore, it is important to carry out a quantitative 
analysis (a) of the structure of this imperfect periodic bending in the bulk, and (b) of 
its influence on the spectrum of undulator radiation. Both of these problems constitute 
the subject of the present paper. 

The paper is organized as follows. In section [2] we briefly outline the basic formulae 
and definitions which refer to a perfect crystalline undulator. Section [3] describes, 
in general terms, the modifications to be introduced to the formalism due to the 
imperfectness of the crystalline undulator. In section H] we present the formalism and 
carry out numerical analysis of the periodic deformations in the bulk caused by a regular 
stress applied to the crystal surface. Finally, in section [5] we carry out quantitative 
analysis of the differences in the emission spectra formed in a perfect undulator and in 
the undulator created by means of periodic surface deformations. 

2. Spectral-angular distribution of the radiation from a perfect crystalline 
undulator 

The spectral distribution of the energy E of radiation emitted by an ultra-relativistic 
(v ~ c) positron in a perfect crystalline undulator can be written in the following form 

d 3 E 

kA ^r. = S(u,6,tp)V N (ri,K d , K a ) , (1) 
ndu d\ I 

where 9 <C 1 and us> are the emission angles with respect to the undulator axis chosen 
along z, dQ = OdOdtp is the emission solid angle. The factor T>n(t], k^, ^a) is explained 
further in this section. The function S(u, 8, tp) is given by 
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Here a ~ 1/137 is the fine structure constant, 7 = e/mc 2 is the relativistic Lorentz 
factor with m standing for the positron mass and e for its energy, ujq = 2nc/\ and p is 
the undulator parameter, 

p = 2vr 7 ^. (4) 
The parameter 77 is defined as follows 

A peculiar feature of the undulator radiation is that for each value of the emission 
angle 9 the spectral distribution consists of a set of narrow and equally spaced peaks 
(harmonics) the frequencies ujj one defines letting 77 = j — 1, 2, 3 . . .: 

The widths Auij of the peaks satisfy the inequality Aujj <C Uj, so that all peaks are well 
separated. 

In an ideal undulator (i.e., in which positrons and photons propagate in vacuumjj] 
the peak intensity is proportional to the squared number of periods. Formally, it 

follows from the fact that d 3 E is proportional to Dn{t]) = (sm.(Niiri)/ sin(7r?7)j which 

behaves as N 2 for integer rj (e.g., [7]). This factor reflects the constructive interference 
of radiation emitted from each of the undulator periods. Consequently, in an ideal 
undulator one can increase unrestrictedly the radiated intensity by increasing of the 
undulator length L = NX. 

The situation is different for a crystalline undulator, where the number of channeling 
particles and the number of photons which can emerge from the crystal decrease with the 
growth of L. In Refs. [T5l HE] the quantitative study of the influence of the dechanneling 
and the photon attenuation on the spectral-angular distribution was presented. The 
main result of this analysis is that the peak intensity is no longer proportional to N 2 . It 
was shown, that in a crystalline undulator the factor D^{t]) must be substituted with 
T>n(jj, Kd, K a ), which depends not only on N and 77 but also on the ratios n d = L/L d and 
n a = L/ L a . Here L d stands for the dechanneling length which is the mean penetration 
distance covered by a channeling particle. The quantity L a , called the attenuation 
length, defines the scale within which the intensity of a photon flux propagating through 
a crystal decreases by a factor of e due to the processes of absorption and scattering. 
A convenient formula for V N (r], k a , /t a ), which enters the right-hand side of ([!]), is as 
follows [IE] : 

4N 2 Ke" Kd 2/t d -fi; a (^ + 40 2 ) e" 

TsNVn, K d, K a ) 



nl + 1QN 2 sin 2 7r(?7 - j v ) [n a - « d « a - K d (2n d - K a ) 2 + A<p 2 
2Kd 2^sm0-(2^ -« )cos^ o _ (2Kd+Ka)/2 l 



(7) 



I The term 'ideal undulator' must not be mixed up with the term 'perfect undulator', which is used 
throughout the paper and stands for the crystalline undulator with a perfect harmonic pattern of 
periodically bent channels. 
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where = 2n{rj — j v )N (with j v being the closest positive integer to rf). 

Despite a cumbersome form of the right-hand side of ([7]) its main features can be 
easily understood. The most important is that, as in the case of an ideal undulator (to 
which T> N (r], K d , K a ) reduces in the limit L d = L a = oo) the main maxima of T> N (rj, K d , K a ) 
correspond to rj integers, and, therefore, the harmonic frequencies are still defined by 
([5]) with 7] = j. For finite Ld and L a , the maximum value of T>]\f(r], k^, K a ) is smaller 
than N 2 whereas the width of the peak is larger than that in the corresponding ideal 
undulator. 

The attenuation length depends on the photon energy, L a = L & (u), and can be 
calculated as the inverse mass attenuation coefficient which are tabulated for all elements 
and for a wide range of photon frequencies [T71 [18] . 

The dechanneling effect stands for a gradual increase in the transverse energy of a 
channeled particle due to inelastic collisions with the crystal constituents [6]. At some 
point the particle gains a transverse energy higher than the planar potential barrier and 
leaves the channel. In a straight crystal dechanneling length Ld of a positron depends 
on the crystal and on the energy of the projectile. In a bent crystal the potential 
barrier changes due to the centrifugal force, and, as a result, Ld becomes dependent 
on the curvature of the channel. A stable channeling of a projectile in a periodically 
bent channel occurs if the maximal centrifugal force in the channel L c f is less than 
the maximal interplanar force £/ max , i.e. C = F c f/U^ ax < 1. For an ultra-relativistic 
particle F c f « e/i? min , where -R m i n is the minimum curvature radius of the bent channel. 
In a perfect crystalline undulator R m in = A 2 /47r 2 a, therefore, the condition for a stable 
channeling reads [2]: 

w max 1 

Using the diffusion model [19] one can be demonstrated that dechanneling length in a 
periodically bent crystal becomes dependent on the parameter C. In the case of ultra- 
relativistic positrons the expression for Ld = Ld(C) can be written as follows [U [20] : 

Ld(C) = (l-CfL d (0), (9) 

where Ld(0) is the dechanneling length a straight channel. This quantity can be 
estimated as Ld(0) = (256/97r 2 )(a,TF dj/ro A c ) 0, H9] . where r and <2tf are the 
classical radius of the electron and the Thomas-Fermi radius of the crystal atom, and 
A c = \\i{^j2rf mc 2 / 1) — 23/24 (/ stands for an average ionization potential of the atom). 

Apart from the dechanneling length, the bending parameter C also defines another 
important quantity, the acceptance of the crystal. The acceptance represents by itself 
the fraction of the particles trapped in the channeling mode (or, in other words, it is 
equal to the ratio of the phase volume corresponding to the channeling regime to the 
phase volume of the beam particles at the entrance of the crystal). In comparison with 
a linear crystal, the acceptance of the bent crystal gains additional decrease due to the 
presence of the centrifugal force in the channels [19]. Generalizing the result of Ref. 
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[2T] . one writes the following expression for the acceptance A(C) of the periodically 
bent channel and for the case of a parallel beam: 

A(C) = (1-C)A(0), (10) 

where .4.(0) is the acceptance of the linear crystal and C < 1. If the centrifugal force 
exceeds the interplanar force (i.e., C > 1) then one uses the identity A(C) = 0, which 
means that the channel does capture the particles. 

In connection with the radiation from a crystalline undulator, the acceptance comes 
into play when one calculates the energy emitted per particle of the bunch. To obtain 
this quantity in a perfect undulator, one multiplies the right-hand side of ([1]) by the 
acceptance A{C). 

The formulae (fTT)- (fT0l) allow one to carry out numerical analysis of the spectral- 
angular distribution formed in a perfect crystalline undulator. In a recent paper [22] 
the results of such analysis were reported for 0.6 and 10 GeV positrons channeling 
through periodically bent Si and Si]_ x Ge x crystals. 

3. Emission from the undulator with a varied bending amplitude 

In this section we analyze the modifications in the formalism of spectral-angular 
distribution of the radiation from a crystalline undulator which appear due to the 
imperfectness of periodic bending of the crystallographic planes. 

For the sake of clarity let us discuss the geometrical conventions which will be 
used below. Figures [1] and [2] help a reader to follow the text. We assume that the 
crystal has the form of a rectangular box. Its length, L, width, I, and thickness, h, 
are measured along the z, x and y directions, respectively. We chose the frame in 
which the values ±h/2 denote the ^/-coordinates of the upper and lower surfaces of the 
crystal. Hence, it is assumed that y = labels the central (xy)-plane of the crystal 
(the midplane). Initially non-deformed crystallographic planes are perpendicular to the 
y axis. We assume that periodic deformation of the crystalline structure occurs only 
in the (yz) -plane, so that there is no deformation in the x-direction. In the deformed 
crystal the bunch of channeling particles propagates in the (yz)-pla.ne along z direction. 
With o y we denote the bunch size along the y direction. The width of the bunch (i.e., 
its size in x) is not used in our model. 

Suppose that bending amplitude is not constant over the crystal thickness but 
changes according to a law a = a(y), as illustrated by figured! The particular form 
of this dependence is discussed in Section H] where we analyze the deformation of the 
crystal interior due to the periodic stress applied to its surfaces. At the moment, the 
only restrictions implied on the dependence a(y) are: (a) the change in amplitude on 
the scale of interplanar spacing d is negligible; (b) a strong inequality a(y) <C A is valid 
for all y. 

Particles in the bunch are randomly distributed along the y-axis. Therefore, we 
may assume that at the entrance a particle can be captured in any channel located 
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within the interval y = [— y max , y max ] > where y max stands for the smallest from h/2 and 
cr y /2. Being captured in the channel at some y point, the particle undulates and emits 
the radiation corresponding to the undulator with the amplitude a(y). For a fixed value 
of the period A, the amplitude defines two important quantities which are the undulator 
parameter p and the bending parameter C. The latter, in turn, defines the dechanneling 
length and the acceptance (see eqs. Qj and ffTUl) ). The influence of p and L d on the 
spectral-angular distribution of radiation is discussed in detail further in this section. 
Here we note that the acceptance A(C) determines the probability of a particle to be 
captured in the channeling more. If C = C(y), this probability depends on the entrance 
coordinate y. Therefore, to obtain the distribution (per particle) of radiation formed 
in the crystal with varied amplitude of bending one should (a) multiply eq. ([1]) by the 
acceptance (fTOl) . and (b) carry out the averaging over the interval y = [— y max , y ma x]- 
This leads to the formula: 



Here d 3 E(y)/hdu dfl stands for the the distribution §1§ obtained for the amplitude 
a(y). The right-hand side of this equation represents the spectral-angular distribution 
per particle averaged over the width of the crystal (or of the bunch if a y < hj^. If the 
amplitude does not change within the interval of integration, the formula (llip reduces 
to eq. ([T]) multiplied by the acceptance corresponding to the fixed C-value. 

Let us discuss in more detail the influence of the dependence a(y) on the 
characteristics of the crystalline undulator radiation. 

It is clear that for the y-varying amplitude the terms, dependent on p in the 
function S(u,9,ip) (see ([2])), change with y since p oc a(y) (see ((4])). Apart from this 
influence, the proportionality can lead to the variation of the harmonics frequencies 
Uj and, consequently, to the loss of the monochromaticity of the radiation. Indeed, 
let a m i n and a max denote the minimum and maximum amplitudes within the interval 
[—Vmax, 2/max]- The corresponding extremum values of the undulator parameter, p min and 
j9 max , having been used in (ED, produce the lower, u^ mm \ and the upper, uj^ max \ bounds 
on Uj. In the limit p max <C 1 there is a weak variation of the harmonics frequencies, 
^jmin) ^ w jmax)^ However, in the case p max > 1 the change in the undulator parameter 

leads to the emission within the band Aoo = u; J - max ' ) — u; J - max ' ) which can greatly exceed not 
only the peak width Autj but also the interval between the neighbouring harmonics. In 
the latter case the monochromaticity of the radiation will be smeared out if one carries 
out the averaging procedure (TTTT) . 

However, the change in the undulator parameter and the harmonics frequencies 
is not the only impact caused by the amplitude variation. As mentioned above, in 
a crystalline undulator the peak intensity is defined by the maximum value of the 
factor V N (r], n d , n a ), which depends on the ratios n d = L/L d and n a = L/L a (see 

§ The contribution to the integral comes only from the regions where C(y) < 1, otherwise A(C) = 0, 




(11) 



see eq. (fTOj) . 



8 



([7])). Typically, in the crystalline undulators based on £ = 0.6 . . . 10 GeV positrons 
channeling in crystals (so far Si monocrystals and Si-Ge mixtures were used [TU [23] ) 
the energy of emitted photons lies within the range 10 2 . . . 10 3 keV. For these energies of 
positrons and photons the attenuation length L a = L a (u) is on the level of several cm 
[TT] and by far exceeds the positron dechanneling lengths (lying within Ld ~ 0.03 . . . 0.7 
cm in straight channels [22j ) . Therefore, assuming K a = in ([7j), one finds that the peak 
value of the factor V N (r], K d , K a ) (i.e., at r\ = j) in an undulator with fixed number of 
periods, N = L/X, is given by 

lim V N (j,K d ,K a ) = 2N 2 l ' {1+ 2 Kd)e - (12) 

The dechanneling length in a channel, which is periodically bent with the amplitude a, 
is defined by equations ([8]) and (Q. Hence, for a = a(y) the right-hand side of ffT2]) (and, 
generally, of ([7]) as well) becomes dependent on y since Kd = L/Ld(C) oc (1 — C)~ 2 with 
C oc a(y). In turn, the variation of the factor D^(rj, k^, K&) can strongly influence the 
averaged spectral-angular distribution (lllj) . 

4. Periodic deformations in bulk 

In this section we present a formalism which allows one to carry out a quantitative 
analysis of the parameters of periodic bending in the bulk of a crystal. We consider the 
case when the bendings in the bulk are due to the periodic deformations on the crystal 
surfaces. Such a situation is illustrated by Figure [2j two parallel opposite surfaces are 
deformed periodically by means of identical sets of parallel grooves applied to each of 
the surfaces (the two sets are shifted by half-period, A/2). 

The deformation of this type can be achieved either by mechanical scratching of 
the crystal surface [T2J [13] or by means of a more accurate laser-ablation method |14j . 
Another possible method is in a deposition of Si 3 N 4 layers onto a Si wafer [13]. 

In either case the crystallographic planes in the bulk become bent periodically 
although the shape y = y(z) of bent planes does not follow an ideal harmonic form 
y = a cos(2tt z / X) . The main deviations are: (a) the amplitude of bending depends 
on the distance from the surface, and (b) in general case, higher subharmonics (i.e. 
the Fourier components of y(z) with smaller periods, X n = X/n, where n = 2,3...) 
contribute noticeably into the formation of the periodic shape. 

However, as we demonstrate below, it is possible to establish the ranges of 
parameters (these include A, the thickness h and elastic constants of the crystal) within 
which the deviations of the resulting periodic shape from the ideal form do not affect 
the spectral-angular distribution of the undulator radiation. 

4-1. The equations of equilibrium with periodic boundary conditions 

To start with, let us formulate the approximations which will be used in our formalism. 
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Figure 2. Schematic representation of a crystal with periodic surface deformations 
(the sets of regularly spaced grooves parallel to the x direction): A stands for the 
period of deformations and 5 denote the width of a groove. The set on the lower 
surface is shifted by A/2 (along the z-axis) with respect to that on the upper surface. 
The surface stress gives rise to the periodic bending of crystallographic planes in the 
bulk of crystal. 



Firstly, we assume that the width of a crystal I (i.e. its size in the x-direction, see 
figure [2]) is much larger than the thickness h: I 3> h. As a result, one can disregard 
the deformations in the x-direction and consider the deformations of a solid body 
in the (yz)-plane only. Additionally, we assume that the period A is incomparably 
smaller than the length L of the crystal in the z direction. Then, taking into account 
the periodicity of the surface deformations, one represents the displacement vector 
u(y, z) = (0,u y (y, z),u z (y, z)) (i.e. the vector which characterizes the change in the 
position vector of a point in the body due to the deformation) in the form of Fourier 
series: 

u(y,z)= s n(y)e mkz , (13) 

71= — OO 

where k = 2tt/\. The vectors s n (y) (with n standing for an integer) are to be defined by 
solving the equation of equilibrium with proper boundary conditions. In what follows 
we adopt that the y coordinate is measured from the crystal midline, and thus y = —h/2 
corresponds to the lower surface, and y = h/2 - to the upper one. 

Secondly, we will consider the limit of small deformations only. In this case the 
strain tensor Uij is defined as follows (see, e.g., |24j ) : 

--K£ + £)' (l4) 

with i = x,y and j = x, y. 

The third approximation concerns the stress tensor, o~ij (we remind, that physically 
<jjj stands for a stress on the i-th plane along the j-th direction). In an isotropic media 
the components of the stress tensor can be related to Uij by means of the two elastic 
constants: the Young's modulus E and the Poisson's ratio v. In the case of a planar 
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deformation the relationship is as follows (see, e.g., §5 in 
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In anisotropic media (e.g., in a crystal) E and v depend on the directions of the applied 
stress and of the deformation. For example, depending on a crystallographic direction 
the Poisson's ratio and the Young's modulus in a Si crystal varies within the intervals 
0.048 < v < 0.403 and 130 < E < 170 GPa [25]. However, to simplify the analysis one 
can chose some average values. In our numerical analysis we use v = 0.28 and E = 150 
GPa which are close to the values used in modeling various deformation processes in 
silicon [13], [26]. To check the sensitivity of the results to the choice of the Poisson's ratio 
we also carried out the calculations using the extreme values of v in silicon. It turned 
out that nearly an order of magnitude change in v does not noticeably affect the results 
which are presented below in the paper. 

The components of the stress tensor satisfy the following equations of equilibrium: 
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dy dz dy dz 

Using (fT3l - (fTB1) in (TIB"]) one derives the system of coupled equations for the functions 



(y) and s nz (y): 
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To find the unique solution of this system one has to impose the boundary conditions. 
These can be formulated as follows. 

Each trench acts as a source of a normal ('J_') and a tangential, or a shear (||) 
tension which are characterized by the average pressures P 1 - and P" applied to the 
crystal surface in vicinity of the trench. In the case when the period A the width 5 (see 
figure [2]), the trenches are equivalent to the sets of concentrated normal and shear forces 
applied along the equally spaced lines on the upper and lower surfaces [jj. Within this 
model the pressures P 1 - and P" can be related to the components of the stress tensor 
calculated at the upper (y = h/2) and the lower (y = —h/2) surfaces. 

Let us first formulate the boundary conditions due to the normal tension. The 
pressure P ± is applied inward the crystal along straight lines, parallel to the x axis, 
passing through the equally-spaced points in the z direction. Therefore, recalling that 
J2j a ij n j represents the 2-th component of the force per unit area (with rij standing 

|| The condition A ^> 6 is well fulfilled in the crystalline undulators which have been manufactured. 
Typical values arc: A = 100 . . . 500 /im and 8 ~ 10 /im [HI EDS [HI [23]. 
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for the j-th component of the outward-pointing normal) one derives the boundary 
conditions: 



a 



mi 



a 



uu 



y=-h/2 



y=h/2 



-XP- 



+00 

E 

n=— 00 
+00 



S(z — An), 



-AP X s(z-\(n+l/2) 



(18) 



a 



11 ~ 



0. 



y=±h/2 

Here the coefficient A ensures that the period- averaged pressure, equals to P -1 . The 
arguments of the delta functions fix the z coordinates of the lines parallel to the y axis. 
As mentioned, the sets of trenches on the upper and lower surfaces are shifted by A/2, 
and this explains the difference of the arguments of the delta functions for y — h/2 and 

y = -h/2. 

In the case of shear stress, supposing that the trench profile is symmetric, one 
notices that the tangential forces created at the opposite edges of a trench are equal in 
magnitude but are antiparallel. Therefore, the component function of z must 

change the sign when crossing a trench. It means that in the limit of infinitesimal width 
5^0 this component becomes proportional not to the delta functions as in (TT8]) but 
to their derivatives. As a result, one writes the boundary conditions in the form: 
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y=±h/2 

where F» = A 2 P" /2% stands for the tangential force associated with the period-averaged 
pressure P". 

With the help of (fT3l)-(fT5l) the equations ( fTHl) and ( fT9|) can be re-written in terms 
of the functions s ny (y) and s nz {y) taken at y — ±h/2. The obtained formulae suffice to 
determine completely the solutions s^f(y) (j = y, z) for the normal or the shear stress. 
To obtain the solution in the case when both types of stress act simultaneously one 
constructs the corresponding linear combinations of the functions s^{y) and s^^y). 

Using these functions further in (|T3l) one determines the displacement vector u(y, z). 
The y component of this vector is of a special interest in connection of the crystalline 
undulator problem since it determines the profile of the periodically bent channel in the 
bulk. 



4-2. Displacement u y (y, z) in the cases of normal and shear stresses 

Let us first analyze the periodic deformation in the bulk due to the normal and the 
shear stresses separately. 
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Resolving the system (j!7p with the boundary conditions (1181) (or (1191) ) one finds the 
functions s^ y (y) (or s ny (y)) for all n. Using these in ffl3|) one represents the y-component 
of the displacement vector in the following form: 

oo 

u^(y,z) = cos(nkz + nn). (20) 

n=l 

Here An (y) can be associated with the amplitude of the nth harmonic of teh periodic 
bending, i.e., the one with the period \ n = \/n. A perfect crystalline undulator 
(see figure [U left panel) is characterized only by the term n = 1 whose amplitude 
is independent on y G [— h/2, h/2]. However, if a crystalline undulator is prepared 
by applying concentrated periodic normal or shear stress then: (a) higher amplitude 
harmonics (with n > 1) appear, and (b) the homogeneity of the bending amplitudes is 
lost since they become ^-dependent. To stress these features of the amplitudes one can 
use the formulae which conveniently expresses An'\y) via v4^'"(0) - the first harmonic 
amplitudes is the midplane of the crystal (i.e. in the (x2;)-plane with y — 0). 
In the case of normal stress alone the formula is as follows: 

{Z^~ cosh(nky) — 2C n nky smh(nky) 
* • " =1 < 3 - 5 < 

smh(nky) — 2S n nky cosh(nky) 
A+ ' 

For a shear stress the relationship reads: 

In cosh(nky) — 2S n nky sinh.(nky) 
A- ' 

n 

Jn smh{nky) — 2C n nky cosh(nky) 



Al(y) = 4(0) ^ x { 
if 



n = 2,4,6... 



n = 1, 3, 5 



(22) 



n = 2,4,6 



A+ 

In (I2T1) and (1221) the following notations are used: 

A± = smh(2<) ± 2n(, = kh fc = 2tt 

C n = cosh(nC), S n = sinh(n^), 2 A 



and 

' = 4(1 - u)C n + 2n(S n , ( X\ = 2(1 - 2v)S n + 2r<C n , 

j£ = 4(1 - u)S n + 2n(C n , \ J\ = 2(1 - 2u)C n + 2<«S n . 

The amplitudes ^'"(O) are proportional to the applied stresses and are given by: 



(24) 



where 



4'"(0) = A^^"(VA) (25) 



2(1 + u) Csinh(C) + 2(1 - u) cosh(C) 

~ sinh(2C) - 2C ' 

2(1 + u) (1 - 2v) sinh(C) + C cosh(C) 1 ' 




sinh(2C) - 2C 




Figure 3. The ratio A^(y)/Af(0) (see (J2TJ) and ([25])) versus y/h (for y e [0,h/2]) 
calculated for several values of the Si crystal thickness as indicated by the parameter 
h/X. The solid curves correspond to n odd (which are indicated explicitly in the top 
left graph), the dashed curves - ton even (see the bottom left graph). 



Let us note that due to the half-period relative displacement of the deformations on 
the upper and lower surfaces of the crystal the amplitudes An (y) with odd values of n 
are even functions of y and vice versa (see (1211) and (I22p ). Therefore, it is sufficient to 
analyze the amplitudes in the upper half of the crystal, i.e. for y = [0, h/2]. 

4-3. Numerical results for the bending amplitudes A^(y) and A n (y) 

The behaviour of the ratios A^(y)/A^(0) and A n {y)/A\{ti) as functions of the scaled 
distance y/h from the midplane is illustrated by figures [3] and HI The dependences were 
calculated for several values of the crystal thickness h (measured in the periods A as 
specified by the ratio h/X) and several even and odd values of n. The data refer to a Si 
crystal and the Poisson's ratio was chosen as v = 0.28. 

It is clearly seen from the figures that for either type of stress the inhomogeneity 
of the bending amplitudes along the y direction is much more pronounced for a thick 
crystal (i.e., when h > X) than for a thin one with h < X. 

Indeed, in the limit h <C A the amplitudes corresponding to odd values of n < X/h 




do not vary noticeably over the crystal thickness. Being even functions of y, the 
odd-n amplitudes An'\y) depend quadratically on y in the vicinity of the midplane, 
An'\y) ~ a + bn 2 (y/X) 2 , and have minimum at y — 0. In a thin crystal the quadratic 
term is small for all y provided n satisfies the condition written above. As a result, 
An'^(y) ~ const for these n. The main reason for the even-n amplitudes to vary even 
in the case of a thin crystal is that they are odd functions of y. Therefore, being 
non-zero at the surfaces y = ±h/2 these amplitudes attain zero in the midplane. (A 
non- monotonous behaviour of A\(y) for particular n and h/X, which is most pronounced 
for the n = 4, 6, 8 curves on the top left graph of figure HI just reflects the fact that the 
right-hand side of (1221) is not, generally, a monotonous function of y.) 

Another important feature of a thin crystal is that over the whole thickness the 
amplitude with n = 1 greatly exceeds those with higher n. This dominance is more 
pronounced for a normal stress than for a shear one as reflected by additional factor n~ l 
on the right-hand side of (I2TI) . As a result, the terms with Af'\y) rs ^'"(0) prevail in 
the series from (fT3l) . so that the periodic bending in a thin crystal is of nearly harmonic 
shape: % (y, z) ~y4^'"(0) cos(27rz/A). 

As the crystal thickness increases the variation of the amplitudes over half-width 
evolves dramatically. In the case of a thick crystal this variation reaches orders of 
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magnitudes. Such a behaviour one understands analyzing the right-hand sides of (1211) 
and (1221) . Assuming h ^> A one derives the following expressions for the amplitudes on 
the surface and in the midplane: 

Ai{h/2) = ^ e l A i(0), Al(h/2) =bjA\(0) 7 n = 1, 2, 3, 4 . . . 
A^(0) = e^~ n ^A\(0), Al(0) =ne^~ n ^ A\(0), n = 1,3,5,... 

where b = (1 — v) in the case of a normal stress and b = (1 — 2v)j2 for a shear stress. 
For even n the identity An (0) = is valid. 

Equations from (12 7p demonstrate that for all n the amplitudes decrease 
exponentially with the penetration distance into the crystal. For a fixed odd n the 
decrease rate can be characterized by the ratio An'\h/2)/An'\o) = be n< >/n(, which 
is, basically, independent on the type of applied stress (the difference manifests itself 
only in the pre-factor b). On the other hand, the first line in (|27|) indicates that in the 
limit of a very thick crystal the amplitudes weakly depend on n in a surface layer of the 
width ~ A (this dependence is more pronounced in the case of a normal stress due to 
the additional factor n). These two features of the amplitude suggest that deviation of 
the periodic bending from the harmonic shape is very strong in the outer layers of the 
crystal whereas in the central layer the terms with n > 1 are negligibly small and the 
profile of bending Uy'\y, z) ~ A± (y) cos(2itz/\) is nearly perfect. The deviation of 
the amplitudes A f'\y) from their values at y = is described by the following formulae: 

2z/-Ctanh(C) 



Af%)KAt\0)(l-^^y where 



K 



2(1 - v) + Ctanh(C) : 
1 + 2z/-Ctanh~ 1 (C) 



(2? 



1 - 2z/ + Ctanh~ 1 (C)' 
To conclude this section we note that figures [3] and 0] suggest that the use of 
crystalline undulators manufactured by means of periodic surface deformation can be 
justified in two cases (irrelevantly to the type of the stress applied to crystal). 

Firstly, it is the limit of a thin crystal h < A which ensures (a) the (nearly) constant 
value of the amplitude Ai(y) m A\(0), and, (b) a small contribution of higher-ra terms 
to sum on the right-hand side of ( |20l) . 

Secondly, for h > A it is meaningful to use only the central part of the crystal, i.e., 
where \y\ < X/2tt, as an undulator. In this case it is necessary to use a narrow beam 
(along the y-direction) which is accurately aligned with the crystal midplane. 

4-4- Estimation of the variation of the undulator parameter p and the bending 
parameter C 

Equation (|20|) shows that in general case the profile of bending contains contributions of 
the terms with various n thus deviating from a pure harmonic shape oc cos(fcz). Hence, 
it is meaningful to analyze the influence of this deviation on the undulator parameter 
p and the bending parameter C which, as it was discussed at the end of section [3], are 
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of importance for the calculation of the averaged spectral-angular distribution fill I) of 
radiation. 

4-4- 1 ■ Variation of the undulator parameter p. Generally, the parameter p characterizes 
the mean-square velocity, Vj_, of the periodic transverse motion of a particle moving 
in an undulator (e.g., [2?]). For an ultra-relativistic projectile the relationship is as 
follows: p 2 = 2 r y 2 v 2 L /c 2 . In a perfect undulator this formula leads to equation (J4]). In 
the imperfect undulator discussed here, a particle moves along the trajectory defined by 
the right-hand side of (120]) with z ~ ct. Calculating the transverse velocity du y (y, z)/dt 
and averaging over the period T = A/c, one derives the expression for the undulator 
parameter as a function of y: 

oo 

p\y) = Y.p 2 n(y^ ( 29 ) 

n=l 

where p n (y) = 27r~fnA n (y) / A stands for the partial undulator parameters corresponding 
to the n-th subhar monies of the trajectory. 




Figure 5. Dependences p(y)/p(0) (see (|29p ) versus y/h calculated for the normal 
(left) and shear (right) stress. Different curves corresponds to several values of Si 
crystal thickness as indicated by the parameter h/X. 



In figure [5] we present the dependences of the ratio p{y)/p{0) on y/h calculated for 
several values of the crystal thickness and for the two types of stress, as indicated in the 
caption. As it was discussed in section [3] the dependence of the undulator parameter 
on y leads to a variation in the harmonic frequencies Uj ([6]) over the crystal thickness. 
If this variation becomes large enough (e.g., when the shift in the harmonic position 
becomes comparable with the separation — ujj of the neighbouring harmonics) 
then the spectral-angular distribution will loose the peak-like pattern typical for the 
undulator radiation. From this viewpoint the graphs from figure [5] allow one to estimate 
the influence of the change in the undulator parameter across the crystal cross section 
on the shape of spectral-angular distribution. 
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For both types of stress the dependences p(y) exhibit the common trend: the 
undulator parameter slowly varies with y in the limit of a thin crystal but becomes 
a rapidly (exponentially) increasing function for a thick crystal. This is a direct 
consequence of the definition (129]) and the behaviour of the amplitudes An(y) (see 
figures [3] and HI). Another peculiarity, which one notices comparing the two panels 
in figure [5], is that for the same value of h/X the ratio p"(y)/p"(0) grows faster than 
p ± (y) /p ± {0). Such a feature, which is more pronounced for higher h/X values, can be 
understood if one recalls the differences in the n-dependences of the bending amplitudes 
in the case of normal and shear stresses, and the behaviour of A n (y) as functions of 
y. Indeed, for moderate and large h/X ratios and for both types of the stress, the 
amplitudes A n (y) with n > 1, being rapidly increasing functions, satisfy the relations 
An(0) < At(0) and A n (h/2) < A^h/2) (see fl27jl ). Hence, for y = the sum on the 
right-hand side of ( [29]) is defined, basically, by the term with n = 1, whereas for larger 
y also the terms with higher n contribute noticeably to the sum. On the other hand, it 
follows from (EQ) and ([22]) that A\(y)/A^(y) oc n, resulting; in a similar estimate for the 
ratio of the partial undulator parameters: Pn(y)/Pn(y) °t n - Therefore, in the case of 
the shear stress the undulator parameter varies with y more rapidly than for the normal 
stress. 

Typically, in crystalline undulators prepared by means of surface deformation and 
based one = 0.6 ... 10 GeV positron channeling [HI [23]), the values of parameter p lie 
within the range ~ 0.1 ... 3. The curves in figure [5] allow us to estimate the degree of 
consistency of using the undulators with various p = p(0) and h/X values. 

It was pointed out in section [3] that there will be no dramatic change in the spectral- 
angular distribution of the radiation due to the imperfectness of the undulator provided 
the variation of the undulator parameter does not lead to a noticeable change in the 
position of the peaks located at the frequencies ouj (see ([6])). To analyze this condition 
let us consider, for the sake of clarity, emission in the forward direction (0 = 0). In this 
case the frequency of the j-th harmonics, radiated from the the channel located at the 
distance y from the midplane, is given by 

^ y) = P^H2- (30) 

It is clear that in the limit p(0) <C 1 the position of the peak will be practically unchanged 
for those y where p 2 (y) <C 1. These two inequalities, depending on the absolute value of 
p(0), allow for a wide-range variation of the ratio p(y) /p(0), up to the order of magnitude. 
It follows from figure [5] that for a thin crystal such a situation can be realizes over the 
full thickness h. In a thick crystal (h > X) only the central part ensures the needed 
variation of the undulator parameter. 

In the opposite limit of large undulator parameters, when p 2 (0) ^ 1, the 
stability of the peak position can be achieved only in that part of the crystal where 
p(0)(p(y) — p(0)J < 1. Taking into account rapid variation of the ratio p(y)/p(0) in the 
case of a thick crystal, one deduces that this inequality can be ensured only by using 
very thin crystals with h <^ X. 
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4-4-2. Variation of the bending parameter C . The channeling process in a bent crystal 
takes place if the centrifugal force, e/R, due to the channel bending is less than the 
interplanar force U' m3X [28]. In a perfect crystalline undulator this condition, applied 
to the points of maximum curvature, results in the inequality (jHJ) which relates {/^ax 
with the period and the amplitude of periodic bending [2j. In the case when the 
profile of periodic bending contains a number of subharmonics, e.g., as in (120]) . is 



it more constructive to relate the parameter C to the mean-square curvature R~ 2 
where the averaging is carried out over the period A. Recalling that the curvature 
R^ 1 is proportional to the modulus of the second derivative of the bending profile, 
\d 2 u y (y, z)/dz 2 \, and carrying out the averaging, one derives the following expression 
for average bending parameter as a function of y: 
/ oo \ 1/2 

c&j= [J2 c n(yn > ( 31 ) 



i 71=1 



where 



1 f ATT 2 en 2 A n (y) 



™ = jW (32) 



is the mean-square partial bending parameter corresponding to the n-th subharmonics 
of the trajectory ( 1201) . 



The parameter C(y), being used in defines the dechanneling length L&(C) as 
a function of y. As mentioned in section [3], the variation of L&(C) over the crystal 
thickness may influence destructively the averaged spectral-angular distribution ( TTTT) . 




Figure 6. Ratios C(y)/C(0) versus y/h calculated for several values of the Si crystal 
thickness (as indicated by the parameter h/X). The left panel presents the dependences 
obtained for the normal stress. The right panel - for the shear stress. 



Figure [6] presents the dependences of C(y), scaled by its value at the crystal center, 
versus y/h calculated for several values of the crystal thickness and for the two types of 
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stress, as indicated in the caption. Qualitatively, the behaviour of the C- L >ll(2/)/C- L 'll(0) 
curves is similar to those of the ratio p ± '^(y)/p ± '^(0) discussed above. Nevertheless, 
there is a quantitative difference: for a fixed value of h/X the curves C- L ^(y)/C ± '^(0) 
increase faster with y than the corresponding ratio of the undulator parameters. This 
feature is due to the difference in the n-dependence of the partial terms in (1201) and in 



(T2"Tj) . Indeed, the partial bending parameters behave as y C%(y) oc n 2 A n (y) (see ( 122}) ) 
instead of the proportionality to nA n (y) of the partial undulator parameters. It was 
discussed in connection with figure [5] that for moderate and large h/X values the terms 
with n > 1 provide the increase of p(y) with y. This increase is more pronounced for 
C(y) since its partial terms contain an extra factor n. 

A crystalline undulator can operate only in the regime when C < 1. If otherwise, 
then the centrifugal force will drive the particles out of the channel. More detailed 
analysis [15] indicated that the reasonable range for the bending parameter is 0.01 . . . 0.3. 
It means that if for y = the bending parameter is of the order of 10 -2 , then the variation 
of C(y) within the order of magnitude are acceptable. From figure [6] it follows that for 
a thin crystal such a situation can be realizes over the full thickness. In a thick crystal 
{h > X) only the central part ensures the needed variation of the bending parameter. 

4-5. Calculation of X corresponding to given amplitudes at the crystal center 

The numerical data, discussed above in sections 14.31 and H~4"t represent the ^-dependence 
of the amplitudes, undulator parameters and bending parameters scaled by their values 
at y = 0. The latter, in turn, can be expressed in terms of the amplitudes A^'"(0) 
which can be calculated from (125]) and (1261) . These equations relate the amplitudes to 
the bending period A. 

In connection with a perfect crystalline undulator it was established (for a detailed 
discussion see [4]) that the operation of the undulator should be considered in the large- 
amplitude regime, i.e. when the bending amplitude is much larger than the interplanar 
distance d. In this limit the characteristic frequencies of undulator and channeling 
radiation (see, e.g., Ref. [29]) are well separated. As a result, the channeling radiation 
does not affect the parameters of the undulator radiation, whereas the intensity of 
undulator radiation becomes comparable or higher than that of the channeling one [3]. 

To apply this approach to a crystalline undulator with varied amplitude one assumes 
that the large-amplitude regime is applicable to the amplitudes at the crystal midplane: 
^'"(0) > d. For the convenience of further consideration let us introduce the quantity 



which explicitly measures the amplitude in the units of interplanar separation. 

Other quantities, which enter equations (1251) and (1261) include: 
(a) The crystalline medium dependent parameters, the Poisson's ratio v and the Young's 
modulus E. (As already mentioned, for a Si crystal we use the average values v = 0.28 
and E = 150 GPa.) 
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(b) The crystal thickness h, which enters via the ratio h/X. 

(c) The applied stress, P 1 - or P'L 

To estimate the stress one takes into account that we are interested in elastic 
deformations of the crystalline structure. Therefore, P^'" must not exceed the plastic 
yield strength, Y, which stands for the stress at which material strain changes from 
elastic deformation to the plastic one. For a silicon crystal one can adopt Y = 7 GPa 
f30] . For further use let us introduce the quantity, which stands for the stress measured 
in the units of Y: 

■■ P J " 



< 1. (34) 

Using ( 133]) and ( 1341) one re-writes equation (125]) as follows (to simplify the notations 
the superscripts and/or '||' are omitted): 
a E d 

x = ^Yj(h/xy (35) 

Explicit forms of the functions J 7 (h/X) for the two types of stress are given in (126]) . 

Relationship (I3"5]) allows one to find the values of A and h which ensure, for a given 
crystal (the parameters E, v and d) and for a relative stress (the parameter k), a desired 
value of the relative amplitude a in the center of crystal. 

To illustrate the relationship we present figure [7J where the dependence of A on 
h/X is plotted for several values of k (as indicated) and for two values of the relative 
amplitude, a = 10 and a = 20. The two panels in the figure correspond to different types 
of the applied stress. It is seen that the A -1 and A" curves obtained for the same values 
of k and a look much alike although there is a distinguishable quantitative difference 
(note the log scale of the vertical axis). 



5. Averaged spectra 

In this section we present the results of numerical analysis of the influence of the 
periodic bending imperfectness on the spectral-angular distribution. The calculations 
were performed for two energies, e = 0.6 GeV and e = 5 eV, of a positron channeling 
along periodically bent (110) crystallographic planes in Si (the interplanar distance 
d = 1.92 A, the maximal interplanar force U^ ax = 6.35 GeV/cm). The data presented 
below refer to the emission in the forward direction, (i.e., 6 = 0° with respect to the 
z axis, see figure [2]). In this case the integrals (j3J) can be evaluated analytically and 
the spectral distribution (pQ) from an undulator with fixed parameters is expressed in 
terms of the Appel function [22], [31] . These features, being not too important physically, 
simplify the numerical analysis. 

The parameters of the crystalline undulators used in the calculations were chosen 
as follows. Firstly, for both mentioned positron energies the amplitude A±(0) of the 
n = 1 subharmonic in the crystal center was fixed as a = A±(0)/d = 10, thus satisfying 
the large-amplitude condition (133]) . Secondly, for each e the length L of the crystal 
was chosen to be equal to the dechanneling length in the straight channel (see (JHJ)): 
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Figure 7. Dependence of A on the relative thickness h/X (see (|35|) ) for several values 
of k = P- 1 -'!! jY (as marked). For each k the solid curve corresponds to A±'^ (0)/d— 10, 
the dashed ones - to Af'"(0)/d = 20 (with d = 1.92 A being the distance between 
the (110) planes in Si). Left panel presents the dependences obtained for the normal 
stress, right panel - for the shear stress. The filled circles and squares on the left panel 
mark the values of A and h for which the averaged spectra (fTTj) were calculated (see 
section [5] for the details). 

L = 0.041 cm and L = 0.31 cm for e = 0.6 and e = 5 GeV, correspondingly. The 
bending period is A = 50 jum for e = 0.6 GeV (resulting in N = L/X = 8 undulator 
periods), and A = 150 /xm for e = 5 GeV (with N = 20). Using ( 1301 with j — 1 one 
calculates the energies of the first harmonics emitted in the perfect undulators (i.e., 
with the cited A and a values): h\o\ = 65.5 keV and h\o\ = 1.2 MeV for e = 0.6 and 
e = 5 GeV, correspondingly. Let us note that the mentioned values of positron energies, 
the crystal lengths and the parameters of periodically bent channels are close to those 
discussed recently in connection with the experiments on crystalline undulators [T4"ll2"2]. 

The results of calculations of the spectral distributions are presented in figures [8] 
and El 

Two graphs in figure [8] correspond to two energies e, as indicated in the caption. 
In each graph, the solid curve represents the profile of the first harmonic peak (in the 
forward direction) calculated using (0Q) for the perfect undulators with the parameters a 
and A given above. Other three curves in each graph correspond to the averaged spectra 
calculated for the same values of a and A but for different crystal thicknesses h. These 
spectra were obtained from ffTTl) by setting y max — h/2 (the effective range of integration 
was restricted by the condition C(y) < 1, see the footnote comment below eq. (TTTT)). 
As it was mentioned in section I4.5[ one can vary the crystal thickness together with 
the relative stress k to achieve fixed values of a (see (1351)). The short-dashed, long- 
dashed and chained curves in the figure were obtained for the relative stress k x = 0.001, 
k 1 - = 0.01 and k 1 - = 0.1, respectively. The corresponding values of h one finds from 
figured (left), where the filled circles mark the ratios h/X for e = 0.6 GeV and the filled 
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photon energy (keV) photon energy (MeV) 



Figure 8. Spectral intensity of the undulator radiation emitted in the forward 
direction by a 0.6 GeV (left panel) and a 5 GeV (right panel) positron channeling 
along periodically bent (110) planes in Si. In each graph the solid curve stands for 
the intensity from a perfect undulator with the fixed amplitudes A\(0) = lOd. The 
dashed, long-dashed and chained curves present the averaged spectra obtained for the 
crystals of different thickness h and exposed to different values of the normal stress. 
Further explanations see in the text. 

squares - for e = 5 GeV. 

Comparing different curves in the figure and recalling the discussion presented in 
section [HI one sees the extent to which the imperfectness of the undulator structure 
over the crystal thickness can influence the emission spectrum. The increase in h (or, 
more generally, in the interval over which the averaging in (|TT|) is carried out) leads to a 
more pronounced variation of the amplitudes A n (y) (see figures [3] and H]) . This, in turn, 
trigger strong variations of the effective undulator parameter p and bending parameter 
C, - figures [5] and El The increase in p with y results in the decrease of the first harmonic 
energy (see (|30p ). In the figure this feature is reflected by the (relative) enhancement 
of the photons with the energy smaller than the values of hwi cited above. Thus, in the 
peaks in the averaged spectrum become wider and the width increases with h. Another 
tendency, clearly seen in the figure, is the decrease in the peak intensity with h. This is 
mainly due to the variation of the bending parameter C(y). Indeed, C varies from its 
minimum value C(0) at the center up to the C(h/2) in the surface layer, figure El For a 
given y the value C(y) defines the dechanneling length L^{C{y)) oc (1 — C(y)) 2 , which, 
in turn, influences the peak value of the spectrum via the factor T>n(t], k&, K a ): smaller 
Ld result in smaller peak intensities (see (JTj) and section [3]). Additionally, the channel 
acceptance (flOl) decreases with C(y) increasing. Therefore, the relative contribution of 
the trajectories with larger C(?/)-values to the integral from (ITTj) increases with h leading 
to the decrease in the peak intensity of the averaged spectrum. 

The two features mentioned above are more pronounced for a positron energy e = 5 
GeV than for e = 0.6 GeV (compare right and left graphs in figure [8]). The reason is a 
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s follows. For a 5 GeV positron, the undulator parameter at the center is p(0) = 0.80 
which is noticeably larger than the value p(0) = 0.28 for a e = 0.6 GeV positron. 
Hence, the influence on the position of the first harmonic peak due to the variation of 
the undulator parameter with y is smaller for lower e. As a result, the averaged peaks 
for e = 5 GeV are (relatively) wider than those for e = 0.6 GeV. The more pronounced 
decrease in the intensities for a 5 GeV positron is due to the larger contribution of the 
trajectories with higher C(y)-values. This, in turns, happens because for the same value 
of the relative stress k the relative thickness h/X is higher for a e = 5 GeV positron 
(compare the positions of filled squares and filled circles which correspond to the same 
k values in figure (0)). 

Figure [8] demonstrates that the pattern of spectral-angular distribution can change 
dramatically. The well-separated peak-like structure, typical for the emission spectrum 
from a perfect undulator, can be completely smeared out in the imperfect undulator due 
to the variation of the parameters of periodic bending. The degree to which the peaks are 
destroyed depends on the crystal thickness and on the values of the undulator parameter 
p(0) and the bending parameter C(0) at the center of the crystal. In particular, for 
p(0) > 1 and in the case of a thick crystal (h ^> A) the averaged spectrum has a nearly 
uniform distribution which is typical for the noise rather than for the undulator-type 
radiation. 

However, even in the case of a comparatively thick crystals one can restore the 
coherence of radiation. To achieve this, it is necessary to avoid using the layers of the 
crystal located far off the centerline, i.e. to use not the whole thickness of the crystal 
but its central part of the width lower than A. In the central layer, the amplitudes 
A n (y) and the related quantities do not deviate noticeably from their values at y = 
(see ( 128]) ). and, therefore, the averaging procedure will not radically influence the peak 
profile. 

To illustrate this statement, we present figure [9] which contains the spectral intensity 
in the vicinity of the first harmonic for the perfect undulator (the solid line) along with 
the averaged spectra calculated for different y max < h/2. The data refer to e = 5 GeV 
positron channeling in a Si crystal with h/X = 1.7, the relative normal stress n 1 - = 0.01, 
the parameters a and A as in figure [HI right. The solid curves in figures [9] and [8] (right) are 
identical. The short-dashed curve (the lowest one) in figure [9] stands for the spectrum 
averaged over the whole thickness, i.e. y ma , x = h/2. Other curves corresponds to smaller 
2/max values, as indicated. It is clearly seen that by narrowing the averaging interval the 
profile of the line can be made closer to that for the perfect undulator. 

6. Conclusions 

In our paper we discussed the influence of imperfect structure of a crystalline undulator 
on the spectral distribution of the radiation. We mainly analyzed the undulators in 
which the periodic bending in the bulk appear as a result of regular surface deformations. 
We demonstrated that this method inevitably leads to two main deviations from the 
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photon energy (MeV) 

Figure 9. Spectral intensity of the undulator radiation in the forward direction 
for a 5 GeV positron channeling along periodically bent Si(llO). The undulator 
period, A = 150 /im, and the crystal thickness, h = 1.7A, correspond to the relative 
normal stress n 1 - = 0.01 (these parameters are indicated in figure [7| by a central 
filled square). The the solid curve stands for the intensity from a perfect undulator 
with the fixed amplitude -Ai(O) = 10d Other curves represent the intensities (fTTj) 
averaged over different intervals y = [0,y max ]: the short-dashed curve corresponds to 
2/max = h/2, the long-dashed curve - to j/ m ax = 0.2/i/2, the short-dashed chained curve 
- to 2/ m ax = 0.1/i/2, the long-dashed chained curve - to t/ max = 0.05/i/2. Further 
explanations see in the text. 

perfect harmonic shape a cos(2ttz/ X). Thsese are: (a) the dependence of the bending 
amplitude on the distance y from crystal midplane, and (b) the presence of subharmonics 
with smaller bending periods. As a result, the quantities which characterize the 
crystalline undulator, - the undulator parameter p and the bending parameter C, vary 
over the crystal thickness h. In turn, this leads to the loss of the monochromaticity of 
the radiation formed in the undulator. 

Typical scale, within which the parameters vary noticeably, is equal to the period 
A of the surface deformations. As we demonstrated in the paper, one can choose the 
following two strategies to partly restore the monochromaticity of radiation. Firstly, one 
can use thin crystals, h < A. In this case, the variation of the amplitude over the width 
as well as the contribution of higher subharmonics do not lead to dramatic changes 
in the spectrum. However, this limit corresponds to very thin crystals, if one takes 
into account that the period of surface deformations lies within the range 50 . . . 200 jivcv 

pa mi. 

The second approach prescribes the use of a thick crystal but in combination with a 
narrow (a y <C A) positron beam along the midplane. This limit seems to be achievable by 
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using existing positron beams in the GeV range and with the size (along one direction) 
of several microns (see [18]). 

To minimize the destructive role of the imperfect structure, one can also consider 
alternative schemes of the surface deformations. Namely, if instead of just a periodic 
surface deformation one applies a harmonic surface deformation with a period A, then 
the only imperfectness of the periodic bending in the bulk will be associated with the 
variation of the amplitude A^y) since for the amplitudes of subharmonics with n > 2 
will be identically equal to zero. This, in turn, will result in a much smaller variation 
of the undulator parameter and the bending parameter over the crystal width (see eqs. 
(129j) and (1311) . and figures [5] and [6]) since this variation is to a great extent due to the 
contribution of higher subharmonics. To achieve the harmonic shape of the surface 
deformation one can either place a crystal between two press molds of the harmonic 
profile (shifted by A/2 with respect to each other) or apply a modulated pressure by 
means of two piezoelectric layers. 
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